Sinusoidal Excitations in Two Component Bose-Einstein Condensates in a Trap 
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The non-linear coupled Gross-Pitaevskii equation governing the dynamics of the two component 
Bose-Einstein condensate (TBEC) is shown to admit sinusoidal, propagating wave solutions in quasi 
one dimensional geometry in a trap. The solutions exist for a wide parameter range, which illustrates 
the procedure for coherent control of these modes through temporal modulation of the parameters, 
like scattering length and oscillator frequency. The effects of time dependent coupling and the trap 
variation on the condensate profile are explicated. The TBEC has also been investigated in presence 
of an optical lattice potential, where the superfluid phase is found to exist under general conditions. 
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Much theoretical work has already gone into study- 
ing the ground state solutions of the coupled Gross- 
Pitaevskii (GP) equations describing multi-component 
BECs [1 H !, Q. TBEC has been observed, where the 
two hyperfine levels of 87 Rb [H, [|| act as the two compo- 
nents. In this fortuitous coincidence in the triplet 
and singlet scattering lengths has led to the suppression 
of exoergic spin-exchange collisions, which lead to heat- 
ing and resultant loss of atoms. A number of interesting 
features, like the preservations of the total density pro- 
file and coherence for a characteristically long time, in 
the face of the phase-diffusing couplings to the environ- 
ment and the complex relative motions Q, point to the 
extremely interesting dynamics of the TBEC. TBEC has 
been produced in a system comprising of K and Kb, 
in which sympathetic cooling of Rb atoms was used to 
condense the K atoms Jj|. It has also been observed in 
7 Li - 133 Cs and 87 Rb - 133 Cs systems [19. 

The presence of nonlinearities in BECs make 
them ideal candidates for observation of solitary waves, 
ubiquitous to non- linear media [l^, [H, Q EH] ■ In the 
TBEC, a number of investigations, primarily devoted to 
the study of localized solitons, have been carried out re- 
cently pi, [13, El, E§ H3|. The coincidence of singlet- 
triplet coupling in 87 Kb, leads to the well known Man- 
akov system 12111 in weak coupling quasi-one dimensional 
scenario [22l . |23| . The rich dynamics of solitons in this in- 
tcgrable system has received considerable attention in the 
literature 0, |H, [H, The effect of spatial inhomo- 



geneity, three-dimensional geometry, and dissipation on 
TBEC have been examined. However, the periodic soli- 
tary waves have not received much attention in the liter- 
ature, particularly in the presence of the harmonic trap 
[28| . Periodic sinusoidal excitations are natural in lin- 
ear systems. In nonlinear models periodic cnoidal waves 
can be present. It is worth mentioning that, in non- 
linear resonant atomic media, cnoidal excitations have 
been experimentally generated p9l . |30| , where relaxation 
naturally led to the atomic level population necessary for 
the existence of these nonlinear periodic waves [3l| . 

Here we analyze the solutions of a generic TBEC model 
in a quasi-one dimensional geometry for periodic solu- 
tions. Interestingly, we find exact sinusoidal wave solu- 
tions in this system in the presence of a harmonic trap, 
which do not occur in the single component case. The 
presence of two components leads to these waves, whose 
energy difference are controlled by the cross phase mod- 
ulation (XPM) . In presence of time dependent trap and 
scattering length, these waves can be compressed and ac- 
celerated. This leads to the possibility of their coherent 
control. We then consider this system in an optical lat- 
tice [H, HH, HH , where a superfluid phase is found to exist 
under general conditions. 

In the case of two species condensate with a wave func- 
tion tf>i(x,t) for the species i, the coupled quasi-lD GP 
equation in the presence of an external potential Vj, can 
be written as, 
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The strength of the intra-species interactions is gi and Vj is the chemical potential. We assume the interspecies 
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form cf)(x, t) = a(t)+b(t)x~\c(t)x 2 tm&l(t) = J Q v(t')dt' . 



interaction to be same for both the components: 512 = Here, <j>(x, t) is a density independent phase having the 
521; Vj is the trapping potential. 

In the absence of any potential, the general traveling 
wave solutions of Eq. (fTa)) and (|lb[) have the following 
form: 



The sinusoidal wave, in this case, is a propagating wave 
with the velocity v(t) in the moving condensate. The 
consistency conditions lead to, 



ipl(x,t) = yt7oi[l-(l ) sin (x - ut)]e 



lp2{x,t) = y£7Q2[l-(l )C0S 2 (X - Ut)] 
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where, ooj's are the equilibrium densities of the atoms in 
the condensed phase. The phase velocity is given by, 
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For these solutions to exist, it is found that the inter- 
actions need to satisfy g 2 2 = 9i<?2 and the background 
densities are related by giOoi = g\i<JQ2- The difference 
between cross phase modulation and self phase modula- 
tion leads to a difference in chemical potentials: 



V1-V2 = (.91 - 912)ctoi(1 + 
= (.912 - .92)0-02(1 + 
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For the limiting case u = 0, the above solutions coincide 
with the solutions mentioned in (33j subjected to the zero 
external periodic potential. 

Recently the effect of the longitudinal trap on the con- 
densate and soliton profile has been investigated quite 
intensively (35j . In the general scenario, the scatter- 
ing length, oscillator frequencies can be time dependent, 
in addition to the presence of a phenomenological loss 
term (36l . l37l . [H, [3^]. Below we employ this method 
to the sinusoidal waves in the two component scenario. 
As will be seen later, this can be used for controlling 
the excitations. They may be compressed or accelerated, 
through suitable temporal modulations of various param- 
eters. We consider self-similar solutions in the oscillator 
trap Vj = \M(t)x 2 , for which the ansatz solution is of 
the following form, 



ipj(x,t) = A{t)a 3 [A{t){x - /(i))] e *[x^*)+<^.*)]. 
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where p, = jUj + Xj. Here Vj{t) = fijA 2 (t) (j = 1,2) and 
Xj 's are constant parameters controlling the energy of the 
excitations. The time dependent wave vector b(t) = A(t) 
and c(t) can be determined by the Ricatti type equation 



^ - k -c\t) = Mit). 
at m 
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From current conservation, amounting to solving the 
imaginary part of the coupled GP equations, one gets 
Eq. ([3]), with the consistency conditions: 



Ut) + —c(t)l(t) - —ACt) = A(t)u 
m m 

= fiA cxp f c / t /w ] 

gj(t) = KjA(t) and g 12 (t) = Ki 2 A(t). 
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FIG. 1: Density distribution of the sinusoidal wave of the first 
component for 87 Kb in presence of a harmonic oscillator trap 
with = 0.4, k 2 = 0.1, u = 0.8 and A = 0.5. 



The real part of the coupled GP equations reduces to, 
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Consistency conditions further require fx = fi = fii + and ki<7oi = ^12002 ■ The form of the densities have been 

Ai = u 2 + A 2 and Ai - A 2 = («i - ki 2 )ct i(1 + 7J jp-) = found to retain their earlier forms: 

2 2 „ 

(^12 — ^2)^02(1 + m fi2 ) with the constraint Aq 2 = ^1^2 
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The non-trivial phases are now controlled by the trap: 



Xi = -r- A(t)[x-l(t)} - tan- 1 !— tan[A(t)(x - l(t))}}, 
n n 
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with a corresponding expression for the second compo- 
nent. The superfluid current densities in presence of the 
trap takes the form 
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with a similar expression for the second component. The 
flow density gets modulated by the chirped phase and as 
expected it depends on the oscillator potential. Hence, 
by tuning the trap the current densities can be controlled 
suitably. 

For illustration, we first consider a trap with M(t) = 
a = const., and inter-species interactions n\ = 0.4 and 
k 2 = 0.1. Mass of the 87 Rb atom is m = 1.41 x I0~ 25 kg. 
The equality of the SPM and XPM leads to the same 
background, along with the same chemical potentials for 
the both the components. Fig. (jTJ) shows the traveling 
wave, with a time dependent velocity in the presence of 
the trap. In presence of oscillator, the atoms can be 
accelerated and suitably controlled. 

It needs to be mentioned that, unlike experimentally 
observed localized solitons, sinusoidal solutions have 
infinite extent, which should be excited in a finite sized 
trap. In a single component BEC, periodic solutions, 
existing in the finite condensate, have been experi- 
mentally seen as Faraday waves [47| . which manifest 
when the scattering length is time dependent in a 
periodic manner [4?J . We expect similar behavior for the 
sinusoidal excitations in two component Bosc-Einstein 
condensates, since these are exact solutions. 

Recently, restricted sinusoidal solutions have been 
found for TBEC in an optical lattice [ijj, where the 
form of the optical lattice potential is taken as, V(x) — 
Vo cos 2 x, where Vq is the amplitude of the optical lat- 
tice. The spatial co-ordinate and Vq are scaled in the 



units of wavelength of incident laser light and recoil en- 
ergy respetively. We find that under general conditions 
the following type of solutions exist: 



1>i(x,t) = ^A + Bcos 2 (x)e tXl ^ +Wlt (13) 
and ip 2 (x,t) = ^C + D cos 2 (x)e lX2(a;)+ ^ 2 *, (14) 

with ujj = | + flj and Xjz = -4f (j = 1,2). Here, c^-s 

are the integration constants. Considering the scenario 
of independent chemical potentials for the two species, 
the consistency conditions yield: 

M2.912 - Hi92 - 2V 2 gi2 „ V 2 gi2 - V\gi , K . 
A = ^ , B = — ^ , (15a) 
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with 
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+ Hj. Dispersion only affects the 



super-current through the integration constants: I\ 
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±j — -n~- The condensate phase for the 
first component has the explicit form: Xi( z ) = 



ci tan [ 



1 r ^A+B tan(z) ■ 



/ \JA{A + B). Similar type of ex- 



pression holds true for the second component. The dif- 
ference between the solutions found here, as compared to 
the earlier one obtained in , lies in the integration con- 
stants U. These constants acquire an additional contri- 
bution from the dispersion term in the form of \AB, not 
present in the restricted solutions found earlier. When 
both the components have identical chemical potentials 
(/ii = U2), the parameter values coincide with Ref. [iol ]. 

In summary, the two component BEC is found to sus- 
tain sinusoidal excitations in a trap, which is not possible 
in the single component case. It is shown that appropri- 
ate changes in the trap and scattering length can be used 
to control the BEC profile. The superfluid velocity can 
also be changed by controlling the experimental param- 
eters. We note that difference between the ground state 
energy of the two components can arise because of the 
XPM. The roles of both harmonic and optical trap to- 
gether is an area worthy of future investigation. It may 
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provide additional parameters for controlling the dynam- 
ical phase transitions found in this system [4l|, Ell, EH . 
One can also study the Faraday patterns in this system 



with time dependent scattering length 0, EH . The pres- 
ence of the two components may affect the nature of these 
excitations. 
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